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Abstract 

We study the Coleman- Weinberg phenomenon in a U (N) x U (N) symmetric scalar 
model in 4 dimensions. By comparing the numerical simulation results with the bare 
perturbation calculation in the weak bare coupling region, we demonstrate explicitly 
that a first order transition is induced by loop-effects. We also observed first order 
phase transitions in the strong coupling region. 



1 Email address: shen@budoe.bu.cdu 



It is well known that a system with more than one relevant coupling can have a first order 
phase transition induced by quantum loop fluctuations - a phenomenon first discovered by 
Coleman and Weinberg |TJ, £| . The application of this Coleman- Weinberg mechanism led 
to a lower bound of 7-10 GeV on the Higgs mass in the Standard Model ||, [5| and to the 
suggestion that the QCD finite temperature chiral restoration phase transition might be 
first order for more than two quark flavor 0. It has also been suggested that the Coleman- 
Weinberg mechanism might be relevant for chiral symmetry breaking in technicolor models 
Recently || it was pointed out that the Coleman- Weinberg mechanism might be used 
to place constraints on the top condensate model or strong-coupling extended technicolor 
models. 



Let us consider the action (in Euclidean form) 

i 2 
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d A x\Ur (d^d^) + ^ftr (0V) + Ai (tr0V) 2 + A 2 tr (fitf] , (1) 



for a scalar field 0, where <p is & N x N complex matrix. Such a model can be considered 
as a low energy effective action of a top condensate model or strongly coupled extended 
technicolor model (For simplicity, we ignore the fermions). It has U(N) x U(N) global 
symmetry under transformation: <p — > U(f)V', where U and V are N x N unitary matrices. 
The dynamical scale A of a top condensate model or strongly coupled extended technicolor 
model provides a natural cut-off. At energies much lower than A the U(N) x U(N) symmetry 
is spontaneously broken to U(N). The vacuum expectation value (VEV) v associated with 
this symmetry breaking is supposed to be responsible for generating the masses of the W 
and Z. For the model in Eq. (|l|) to be a good low energy effective theory, there must be a 
large hierarchy between the VEV scale v and the cut-off scale A. 

However, this model is expected to have first order phase transitions because of the 
Coleman- Weinberg mechanism [0, ||. As m? is tuned through the critical point, the VEV 
jumps from zero in the symmetric phase to a finite nonzero value in the broken phase. In 
certain portions of parameter space the transition may become strongly first order such that 
the condition v/A -C 1 can no longer be maintained. Therefore this part of the parameter 
space should be excluded from the low energy effective model Eq. (|I]). It is hoped that this 
in turn will constrain the parameter space of the top condensate model or strong-coupling 
extended technicolor model M. 

For the convenience of later discussions we first reproduce a geometric description of the 
Coleman- Weinberg mechanism developed by Yamagishi JT0], which combines the effective 
potential approach originally used by Coleman and Weinberg and later developments in 
terms of renormalization group (RG) flows and fixed point structure in parameter space 
1,0,1- 

In the parameter region X 2 > 0, spontaneous symmetry breaking happens along the 
diagonal of matrix with a minimum at (j)^ = ±v5ij. The tree level potential in a zero-mass 
theory can be written as 

Ufa) = N(NX 1 + A 2 V 4 , (2) 



where ip = Jtrfaftcf)) /N . Away from ip = quantum fluctuation will contribute. However, 
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according to the solution of the RG equation for U(<p), the loop corrections can be included 
by letting Ai, A2 run with ip while keeping the form of U(ip) in Eq. (0) unchanged (here we 
have accounted the fact that the anomalous dimension, 77, is zero in leading order for this 
model). If we define t = ln((p/n) with /1 a renormalization scale where Ai,A2 are defined, 
the running of Xi(t), \2(t) will be governed by 3 functions 

^ = A(Ai,A a ), *± = p a (\ 1 ,\ 2 ) (3) 

with initial conditions Ai(0) = Ai,A2(0) = A2. A local minimum will develop and the 
Coleman- Weinberg phenomenon will be found if the running of Ai, A2 crosses the "stability 
line" 

4(iVAi + A 2 ) + NPi + p 2 = , (4) 

in a region where 

NX, + A 2 > , 4(NP 1 + 3 2 ) +iVl>tr + > • (5) 



In the region A2 < the symmetry breaking happens along one diagonal element fcj = 
±v5u5ij. In this case we write the tree level potential as 

[%) = (A! + A2V 4 , (6) 

with ip = ^Jtr (ftcj)). Again the loop corrections can be included in the running of Ai,A2- 
The Coleman- Weinberg phenomenon will occur if the running of Ai, A2 crosses the "stability 
line" 

4(A 1 + A 2 )+/3i + /?2 = 0, (7) 

in a region where 
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A 1 + A 2 >0, 4(A + ft)+^ft^>0, (8) 



Using the one-loop perturbative /5-functions 0, |^| 

(^ 2 + 4)A 2 + 4ATA 1 A 2 + 3A 2 , (9) 
6A X A 2 + 2A^A 2 , 

the RG flow of Ai, A2 is plotted in Figure 1 for N = 2. Notice that the arrows in Figure 1 
indicate the direction of decreasing ip (i.e. the infrared limit of the theory). The "stability 
lines" are indicated by the dotted lines. They are entirely inside the regions that satisfy 
the conditions of Eqs. and Eqs. (^). Thus whenever the flow lines cross the "stability 
line" one expects the Coleman- Weinberg phenomenon^. In the language of the RG flow the 

2 In Figure 1, the perturbative "stability lines" terminate in the strong coupling region (solid circles). 
This may merely mean the breaking down of perturbation theory. Indeed, we find in numerical simulation 
a first order phase transition at Ai = —19.0, A2 = 40.0. 



2 dXi 

~dT 
2 dX-2_ 
~dT 
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Coleman- Weinberg phenomenon occurs due to the absence of an infrared stable fixed point. 
The system is driven to a runaway RG trajectory in the infrared limit and the vacuum at 
(p = becomes unstable |3|, [7], |^| . 

The above description is simple and elegant. However, the /3-functions and consequently 
the RG flow in Figure 1 was calculated in the leading order renormalized perturbation theory. 
In principle it is possible that there might exist a nonperturbative infrared stable fixed point 
and the phase transition could become second order inside its attractive domain. For the 
purpose of Ref f|, this is an important question. If one stays around Ai = A2 = the 
couplings run very slowly. It takes many decades to reach the VEV scale from the cut-off 
scale. Then the phase transition becomes practically second order and there can be a large 
hierarchy between A and the VEV. It is only in the region where Ai is small and A2 is 
relatively large that the RG flow becomes fast and a small A/v ratio becomes likely. This 
has been used |§ as an argument to exclude the large A 2 region from the model in Eq. (|l|). 
However, should there be an infrared stable fixed point in the large parameter region, the 
model in Eq. (TJ) would have practically no constraints on its parameter values (besides 
the obvious condition A^Ai + A2 > for the bare action to be stable). Therefore, it was 
suggested || that a numerical simulation is needed to look for first order phase transitions 
in the strong coupling region. 

Interest in numerical studies of the Coleman- Weinberg phenomenon has also been raised 
in a recent work of March- Russell [TT|] . He considered gauge-Higgs systems in 3 dimensions 



which are thought to have first order phase transitions due to the Coleman- Weinberg mech- 
anism. He argued that the conventional picture based on the d = 4 — e expansion result 
is unreliable for e = 1 and the phase transition for such systems in 3 dimensions could be 
second ordeif]. The model in Eq. (fl|) can serve as a testing ground for this conjecture []. 

In 4 dimensions there is no problem associated with the e-expansion and the Coleman- 
Weinberg phenomenon is expected to exist at least in the weak coupling region. However, 
clear numerical evidence for such phenomenon has not been obtained so far. In the follow- 
ing, we report our first results in the numerical investigation of the model in Eq. (|XJ) in 4 
dimensions. 

For the numerical simulations we set N = 2 in Eq. (TI]). The simple Metropolis algorithm 
is used with a uniform step size A0 i:) = A0 which is tuned such that the acceptance rate is 
around 60%. Since there can be no symmetry breaking in a finite volume it is important to 
choose an order parameter invariant under U(N) x U(N). For a lattice with linear size L, 
we define 

4>v = JJ[12M X ) , (10) 



3 Numerical simulations of a spin model (T^] and SU(N), U(N) nonlinear sigma models |l3) in 3 dimensions 
found evidences of first order phase transitions. These works seem to favor the e-expansion predictions. 

4 A Monte-Carlo simulation had been performed before for the model in Eq. (|l|) in 3 dimensions J[4| . 
Unfortunately the result of Ref. |Q was wrong. The point shown to have two coexisting phases in that 
paper (Figure 2) actually is deep in the broken phase. Our investigation of this model in 3 dimensions found 
first order phase transitions and the results will be reported elsewhere. 
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and use (tr(0V)) as the order parameter. This order parameter is a finite number in the 
broken phase and vanishes like 0{1/L A ) in the symmetric phase. 

For this exploratory investigation most simulations were done on relatively small lattices 
(4 4 , 6 4 and 10 4 ). To decide the order of phase transitions we looked for the hysteresis effects 
in the thermocycles: we perform a series of runs back and forth across the critical region. 
Each run uses the last configuration of the previous run as its initial configuration. In case 
of a first order phase transition, supercooling or superheating will produce the hysteresis 
effects. In addition, we looked at the histogram for the distribution of £r(</>'</>). Since 
the histogram distribution is proportional to exp{— L 4 Ul{4>)} , where Ul(4>) is the effective 
potential, one should observe a double peak structure throughout the critical region in a first 
order phase transition. There are more elaborate techniques [|15] to determine the order of 
phase transitions which will not be used here. 

We choose the parameters Ai and A2 such that the bare lattice action is absolutely stable 

A 2 >0, iVAi + A 2 >0. (11) 

According to the RG flow in Figure 1, Ai and A 2 run very slowly in the weak coupling 
region. Qualitatively, a weak coupling bare action fixed to the right of the "stability line" in 
Figure 1 , will need many decades of running in energy scale to get a renormalized theory to 
cross the "stability line". The running becomes fast in the strong A 2 region. On the lattice 
the running is controlled by the correlation length £ with t « ln£, which is limited by the 
lattice size £<L. Therefore, to observe the Coleman- Weinberg phenomenon one either has 
to perform simulations on a large enough lattice 5 or has to move the bare theory very close 
to the "stability line" in order to get a renormalized theory to cross the "stability line" . 

Our numerical results in the weak coupling region are shown in Figures 2 and 3 where 
the bare parameters are chosen to be very close to the "stability line" Eq. ([|). The phase 
transition is clearly first order at Ai = —0.22, A 2 = 0.5. Figure 2 shows thermocycles for the 
4 4 , 6 4 and 10 4 lattices. For the 4 4 lattice each point in the thermocycle has 20000 sweeps as 
warm-up and 60000 sweeps in the measurement. For the 6 4 and 10 4 lattices they are 3000, 
10000 and 400, 1200, respectively. Away from the critical point the finite size effects are 
small. They become larger close to the critical point: the hysteresis effect is barely visible on 
the 4 4 lattice and becomes stronger on larger lattices as shown in Figure 2. Figure 3 shows 
the tunnelings between the broken and symmetric phases on the 4 4 lattice at the "critical 
point" . 

In the weak coupling region one may calculate the order parameter (trffiif))) in the bare 
perturbation expansion and compare with the numerical results. In particular, for A 2 > the 



effective potential Ul{4>) has a global minimum at <pij = ±vdij. If we define tp = ytr((ffl(f))/N, 
the one-loop effective potential is given by 

U(<p) = A^^y? 2 + (iVAx + A 2 )y/) + ^^[p 2 + m 2 + U^X, + \ 2 ) V 



(12) 



3 We believe this might be the physical reason for the phenomenon observed in Ref. pi 
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where d = 4 and p 2 = X^sin 2 ^. In the infinite volume limit the lattice sum should be 
replaced by a momentum space integral J2 P ^o I (f^p- 

At tree-level the effective potential Eq. (12) appears to indicate a second order phase 
transition at m 2 = for arbitrary values of Ai, A2 (provided A^Ai + A2 > 0). This is, 
however, changed at one-loop. For example, if we plot the one-loop effective potential at 
Ai = —0.22, A2 = 0.5 for various m 2 values. In contrast to the tree-level picture where there 
is always a single minimum, we will have two minima at (p — and (p — v > 0. As m 2 is 
tuned from —1.0 to —1.1 the second minimum at <p > changes from a local minimum to a 
global minimum: a picture ordinarily found at a first order phase transition. 

The dotted line in Figure 2 indicates (ir(0'0)) as calculated in the bare perturbation 
theory. The agreement with the numerical result is reasonably well considering that A2 = 0.5 
is not exactly weak coupling with our normalization. Numerically we found that the bare 
perturbation theory works within 20% up to Ai, A2 ~ 1. 

Since very often first order phase transition happens at a finite correlation length where 
the cut-off effects can not be neglected (known as "no continuum limit"), one may wonder 
why the description of the Coleman- Weinberg phenomenon in the renormalized language 
could be qualitatively correct. The answer lies in the fact that the Coleman- Weinberg 
phenomenon happens in the infrared. Given a renormalized theory \m, X 2 r fixed to the left 
of the "stability line" in Figure 1, one can always arrange (at least in the small coupling 
region) the bare theory sufficiently far away from the "stability line" along a trajectory such 
that it takes many decades in energy scale for the bare theory to "run" into the renormalized 
theory. In such arrangement the cut-off effects can be made small (the phase transition 
becomes weakly first order) and the use of the renormalized theory language becomes valid. 

Of course, the cut-off effects can not be neglected in a parameter region where the phase 
transition is strongly first order. As pointed out in Ref. ||, such region must be excluded 
from the continuum theory because the renormalized theory is no longer meaningful^. 

In the strong coupling region we also found first order phase transitions. For example, in 
Figure 4 we show our simulation results at Ai = 0, A2 = 10 on the 10 4 lattice. The lower and 
upper curve correspond to Monte Carlo runs with disordered and ordered initial conditions. 
The fact that we are able to observe a first order phase transition here on a relatively small 
lattice conforms to the qualitative feature of the RG: although the point Ai = 0, A 2 = 10 
appears to be far away from the "stability line" in Figure 1, the running is fast due to the 
strong A 2 coupling. 

In conclusion, by comparing the numerical simulation results with the one-loop bare 
perturbation calculation in the weak coupling region, we have shown explicitly that in d = 4 



6 This is similar in spirit to the argument for the triviality bound on the Higgs mass fllq] . To actually map 
out the excluded region in the (Ai, A2) parameter space requires detailed work and is under investigation. 
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quantum fluctuations can induce a first order phase transition-commonly known as the 
Coleman- Weinberg phenomenon ]1] . We also found first order phase transitions in the strong 
A2 coupling region. This might be taken as an indication that there is no infrared stable 
fixed point in the nonperturbative region. However, a systematic search (which is under 
investigation) is needed in order to confirm this point. 

The numerical work to determine the complete phase diagram in both weak and strong 
coupling regions will be reported elsewhere |T7| . 
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Figure Captions 



Figure 1: Solid lines indicate the RG flow in Ai, A2 parameter space. The arrows point to 
the infrared direction. The dotted lines are the "stability lines" calculated according to Eqs. 
(£D and ([!]) using the one-loop /3- functions. They stop at end points (solid circle). 

Figure 2: (trffify) as a function of m 2 at Ai = — 0.22, A2 = 0.5. The open circles indi- 
cate data on the 4 4 lattice and open squares and open triangles are for 6 4 and 10 4 lattice 
respectively. Unless explicitly shown the statistical errors for data points are smaller than 
the size of symbols. The solid lines connect data points from 6 4 and 10 4 lattices to indicate 
the thermocycles. The dotted line gives the bare perturbation prediction. 

Figure 3: Time history shows tunneling between two states at Ai = —0.22, A2 = 0.5, m 2 = 
—0.875. Data is obtained on the 4 4 lattice. Tunneling becomes rare on larger lattices. 

Figure 4: Time history shows two-states signal on the 10 4 lattice at Ai = 0, A2 = 10, m 2 = 
—29.625. The lower and upper curve correspond to runs with disordered and ordered initial 
conditions. 
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